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1.0  INTRODUCTION 


-*  Ionograms  contain  as  basic  ^information  the  virtual  range 
versus  frequency.  The  Digisonde  128PS  [Bibl  and  Reinisch,  1978]  also 
measures  amplitude,  phase,  polarization,  incidence  angle  and  the 
spectral  signature  of  the  echoes.  This  allows  automatic  scaling  of 
the  ionograms  and  conversion  to  vertical  electron  density  profiles. 
This  report  discusses  the  automatic  scaling  of  some  8000  Digisonde 
ionograms  from  the  Goose  Bay  Ionospheric  Observatory  (64. 6N 
geomagnetic)  of  the  Air  Force  Geophysics  Laboratory.  These  ionograms 
display  a  large  variety  of  features:  quiet  and  disturbed  daytime 
recordings,  heavy  spread  F  during  the  night,  the  mid-latitude  trough 
moving  over  the  station,  fast  variation  of  the  ionospheric  parameters 
and  frequent  absorption  events. 

Sections  2,  3  and  4  give  a  brief  description  of  the  data 
base  and  the  automatic  scaling  methods  that  extract  the  required 
ionospheric  parameters  from  the  digital  ionograms.  A  few  examples 
illustrate  the  results  of  the  profile-fitting  method  that  is  used  to 
calculate  the  electron  density  profiles  from  the  autoscaled  echo 
traces.  It  is  the  first  time  that  anybody  has  succeeded  to 
automatically  scale  a  large  number  of  quiet  and  disturbed  ionograms." 
To  establish  a  level  of  confidence  in  the  scaling  algorithm,  the 
critical  frequency  and  height  values  for  more  than  2000  hourly 
ionograms  of  January,  April,  July  and  September  1980  were  compared 
with  the  manually  scaled  values.  Section  5  shows  that  for  90%  of  all 
ionograms  the  error  in  foF2  is  less  than  0.5  MHz;  the  rate  is  95%  for 
a  1  MHz  error  limit.  This  is  an  excellent  result  considering  the 
subauroral  location  of  Goose  Bay.  The  Automatic  Real  Time  Jonogram 
Scaler  with  True  Height  Calculation  (ARTIST)  which  is  based  on  the 
8086  microprocessor  is  briefly  described  in  Section  6. 
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2.0  DATA  BASE 


The  routine  ionograms  at  the  Goose  Bay  Ionospheric 
Observatory  (53. 3N,  60. 5W  geographic,  64. 6N,  12. 1W  geomagnetic) 
measure  and  record  amplitude,  polarization,  incidence  angle  and 
Doppler  frequency  as  a  function  of  frequency  and  range.  The  lower 
part  of  Figure  1  shows  a  quiet  amplitude  ionogram  containing  all 
signals.  Removal  of  the  non-vertical  and  X-polarization  signals 
results  in  the  upper  ionogram  of  Figure  1,  which  is  much  easier  to 
scale  automatically.  The  X-trace  data  are  not  discarded;  they  are 
used  for  the  accurate  determination  of  foF2.  For  bottomside 
ionograms,  the  0-trace  is  generally  better  presented  than  the  X-trace 
and  our  autoscaling  effort  concentrated  therefore  on  the  0-trace. 
This  is  in  contrast  to  topside  ionograms  where  0  and  X-polarization 
echoes  dominate  alternatingly  [Reinisch  and  Huang,  1982]. 

Four  months  of  tape  recorded  ionograms  for  January,  April, 
July  and  September  1980  were  processed  to  test  for  any  seasonal 
dependence  of  the  scaling  accuracy.  Three  vertical  ionograms  were 
recorded  per  hour,  except  in  September  when  6  ionograms  per  hour  were 
stored.  Disregarding  any  of  the  technically  incomplete  ionograms, 
about  8000  ionograms  were  autoscaled.  The  parameters  of  2226  hourly 
ionograms  were  compared  with  the  manually  scaled  values. 
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Figure  1 
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3.0  SCALING  ALGORITHM 


The  ideas  and  procedures  for  the  automatic  scaling  of 
digital  ionograms  have  been  published  in  a  series  of  papers  [Reinisch 
and  Huang,  1982;  Huang  and  Reinisch,  1982;  Reinisch,  Moses  and  Tang, 
1981],  A  paper  on  the  processing  of  bottomside  Digisonde  ionograms 
has  been  submitted  for  publication  to  Radio  Science.  Appendix  A  is  a 
preprint  of  this  paper,  describing  in  detail  the  scaling  algorithms 
and  the  profile-fitting  method  developed  for  the  calculation  of  the 
electron  density  profile.  It  suffices  to  give  here  a  brief  summary. 

Ideal  ionograms,  recorded  under  quiet  conditions  with  a 
relatively  low  level  of  interference,  pose  no  difficulties  for 
automatic  scaling,  yet  they  are  useful  to  illustrate  some  of  the 
scaling  procedures.  (Manufacturers  of  ionosondes  generally  select 
such  quiet  ionograms  to  demonstrate  system  performance  and  the 
capability  of  automatic  scaling.  We  refer  to  Figures  3,  4  and  7  in 
Appendix  A  for  examples  of  disturbed  ionograms  which  are  successfully 
scaled  by  the  ULCAR  algorithm.)  It  is  important  to  emphasize  that 
even  at  mid-latitude  the  percentage  of  unusual  or  disturbed  ionograms 
might  be  as  high  as  30%.  The  basic  concepts  of  the  scaling  procedure 
must  take  into  account  this  fact.  In  general,  the  vertical  0  and  X 
echo  traces  must  be  found  within  spread  F  signals.  Multiples  and 
oblique  echoes  must  be  eliminated  relying  on  the  amplitude,  polariza¬ 
tion  and  incidence  angle  information  contained  in  the  Digisonde  iono¬ 
grams.  This  is  only  possible  by  examining  the  ionogram  in  its  entirety. 

To  find  the  F  trace  the  ionogram  is  surveyed  for  heights 
larger  than  160  km,  and  the  "center  window"  (Figure  1)  with  the 
maximum  signal  energy  is  determined.  A  first  approximate  trace,  the 
baseline,  is  constructed  by  sliding  a  searching  window  toward  higher 
and  lower  frequencies.  Two  hyperbolic  functions 

h'x  *  r  +  a+ET 


h' 


0 


and 


r  +  a+blT+ifHj 


are  fitted  to  the  0  and  X  amplitude  pixels  in  ordtr  to  determine  foF2, 
the  critical  frequency  of  the  F2  layer;  f  and  fH  are  sounder  and 
gyrofrequency,  respectively.  The  block  diagram  in  Figure  2  gives  an 
outline  of  the  F-trace  extraction. 

The  main  difficulties  for  the  E-region  are  the 
identification  of  £,  Es  and  night  E.  To  find  the  normal  E  trace  an 
analytic  function  is  fitted  to  the  amplitude  pixels.  The  function  is 
derived  from  a  parabolic  profile  and  the  three  parameters  of  the 
parabola,  height,  half  width  and  peak  density,  are  determined  such  as 
to  maximize  the  average  signal  amplitude  of  the  ordinary  vertical 
echoes  traced  out  by  the  h ' (f )  function.  Continuous  0  echoes  beyond 
foE  are  identified  as  Es  trace.  To  save  CPU  time,  the  search  for  the 
peak  density  (or  foE)  of  the  parabolic  E-layer  is  limited  to  ±0.3  MHz 
around  a  predicted  mean  value.  For  evening  and  night  times  the  search 
is  extended  to  6  MHz,  an  empirical  upper  limit  for  night  E  critical 
frequencies. 


4.0  SCALED  PARAMETERS  AND  ELECTRON  DENSITY  PROFILES 


Once  the  E  and  F  region  traces  h ' (f ^ )  are  focused,  the 
program  determines  the  following  parameters. 


(1)  foF2 

F2  layer  critical  frequency 

(2)  foFl 

FI  layer  critical  frequency 

(3)  MUF(3000) 

Maximum  usable  frequency  for  D  =  3000  km 

(4)  M(3000) 

MUF  (3000)/foF2 

(5)  fminF 

Minimum  frequency  for  F  echoes 

(6)  fmin 

Minimum  frequency  for  E  or  F  echoes 

(7)  foE 

E  layer  critical  frequency 

(8)  foEs 

Es  layer  critical  frequency 

(9)  h'F 

Minimum  virtual  height  of  F  trace 

(10)  h'F2 

Minimum  virtual  height  of  F2  trace 

(11)  h'E 

Minimum  virtual  height  of  E  trace 

(12)  h'Es 

Minimum  virtual  height  of  Es  trace 

The  scaling  accuracy  of  parameters  (1)  to  (12)  is  statistically 
evaluated  in  Section  5.  The  following  parameters  are  also  determined 
but  no  evaluation  is  given  in  this  report. 

(13)  FQ 

Average  range  spread  of  F  echoes 

(14)  fXI 

Maximum  frequency  of  F  trace;  fXI-fXF2=F  region 
frequency  spread 

(15)  EO 

Average  range  spread  of  E  echoes 

(16)  fminE 

Minimum  frequency  of  E  echoes 

(17)  A 

Echo  amplitudes,  normalized  to  100  km  reflection 
height,  for  each  MHz 

(18)  Vz 

Vertical  velocity  deduced  from  doppler 

measurements 
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All  these  parameters  are  available  as  program  output  together  with  the 
h ' ( f ^ )  array. 

Having  h ' ( f ^ )  available  it  is  of  course  possible  to 
calculate  the  electron  density  profile  [for  example,  Doupnik  and 
Schmerl ing,  1965].  Since  the  autoscaling  method  produces  some 
occasional  wild  points  which  are  likely  to  adversely  effect  the 
conventional  lamination  procedure  we  developed  the  profile-fitting 
technique  for  the  calculation  of  the  F-region  profile.  The  F-profile 
is  represented  by  a  single  analytical  function  consisting  of  a 
modified  sum  of  shifted  Chebyshev  polynomials  [Huang  and  Reinisch, 
1982,  p.  838,  eq.  6].  The  idea  of  polynomial  fitting  had  been 
suggested  earlier  [Titheridge,  1967].  A  parabolic  approximation  of 
the  E-region  profile  is  automatically  obtained  in  the  E  layer  scaling 
routine.  E  and  F  profiles  are  joining  smoothly  allowing  for  a 
parabolic  valley  (see  Figure  9  in  Appendix  A)  between  the  E  and  F 
region.  Comparisons  between  manual  and  autoprofiles  are  given  in 
Figures  10,  11  and  12  in  Appendix  A. 
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5.0  PERFORMANCE  EVALUATION  OF  THE  IONOGRAM 


SCALING  ALGORITHM 


It  is  always  easy  to  automatically  scale  a  few  selected 
ionograms  by  tailoring  the  scaling  algorithms  to  these  specific 
ionograms.  The  real  test  is  the  application  of  the  algorithm  to  a 
large  number  of  ionograms  covering  day,  night  and  twilight  as  well  as 
quiet  and  disturbed  conditions.  To  test  the  performance  of  the  BISA 
(Bottomside  Ionogram  Scaling  Algorithm)  program,  more  than  8000 
ionograms  for  January,  April,  July  and  September  1980  from  Goose  Bay, 
Labrador,  were  processed.  This  data  base  is  representative  of  one 
year's  data  covering  all  seasons  and  all  types  of  ionospheric 
conditions.  The  parameters  of  the  approximately  2200  hourly  ionograms 
were  compared  with  the  manually  scaled  values. 

The  difficulty  of  the  scaling  task  is  illustrated  by  the 

sequence  of  ionograms  in  Figure  3  recorded  in  Goose  Bay  on  7  January 

1980  between  1819  and  2039  AST.  It  is  the  time  when  the  F-region 
trough  moves  over  the  station  causing  spread  F  and  oblique  echoes. 
The  preliminary  results  of  the  autoscaling,  indicated  by  the  small 
filled  circles,  are  in  good  agreement  with  manual  interpretation. 

(The  final  h'(f)  traces  are  obtained  by  lowering  the  preliminary 
traces  to  best  fit  the  leading  edge  of  the  pulse.)  The  amplitude 

ionograms  in  Figure  3  contain  the  0  and  X-echoes  as  well  as  oblique 
returns.  If  reflections  from  the  trough  walls  or  other 
irregularities,  like  at  1959  AST,  produce  several  retardation  cusps  in 
the  ionogram,  the  program  selects  the  vertical  0-echo  sequence  with 
the  strongest  amplitudes.  A  manual  evaluator  may  scale  foF2  at  1959 
AST  as  2.7  or  3.2  MHz,  while  the  program  scaled  4.3  MHz  and  jumped  to 
the  lower  value  at  2019  AST.  These  differences  in  "interpretation" 
are  the  main  contributor  to  errors  in  foF2  above  0.5  MHz. 


5.1  Comparing  Manual  with  Autoscaling  Results 

The  parameters  obtained  from  the  BISA  program  are  compared 
to  the  corresponding  manually  scaled  parameters  to  test  the  accuracy 
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of  autoscaling.  For  2226  ionograms  the  parameters  foF2,  M(3000), 
MUF(3000),  foFl ,  h'F,  h'F2,  fminF,  fmin,  foE,  foEs,  h'E  and  h'Es  were 
investigated  for  their  accuracy.  The  BISA  program  separates  the 
E-region  parameters  into  day  and  night  groups  to  independently  assess 
the  accuracy  of  night  E  scaling.  Ionograms  with  technical  errors  were 
removed  from  the  data  base. 


5.1.1  foF2 

The  critical  frequency  of  the  F2  layer  is  perhaps  the  most 
important  ionospheric  parameter.  The  minimum  accuracy  requirement  for 
a  high  latitude  station  like  Goose  Bay  was  set  to  ±1  MHz  for  80%  of 
all  ionograms.  Figures  4  to  7  show  the  error  distributions  for  the 
selected  four  months.  It  can  be  seen  that  the  minimum  requirement  is 
by  far  exceeded.  The  1  MHz  error  limit  test  is  passed  by  92%  of  the 
ionograms  in  January,  96%  in  April,  97%  in  July  and  in  September. 
Indeed,  90%  of  all  ionograms  satisfy  the  0.5  MHz  error  limit.  The 
curves  on  the  left  side  of  the  figures  show  the  symmetry  of  the  error 
distribution  which  is  more  or  less  Gaussian  in  shape.  The  error 
statistics  for  the  ionograms  without  spread  F,  e.e.,  daytime 
ionograms,  were  established  separately.  It  is  satisfying  to  see  that 
the  error  curves  for  all  ionograms  (dashed  curves)  are  only  about  10% 
lower  than  those  for  non-spread  ionograms  (solid  curves). 


5.1.2  MUF(3000) 

Figures  8  to  11  display  the  error  distribution  function  of 
MUF(3000)  where  the  error  is  defined  as  the  percentage  difference 
between  the  manual  and  autoscaled  value  based  on  the  manual  reading. 

In  July  and  September,  over  90%  of  all  ionograms  have  less  than  10% 

errors,  in  April  88%  and  in  January  80%.  Averaged  over  the  four 

months,  88%  of  all  ionograms  are  scaled  with  a  MUF(3000)  value  of  less 
than  10%.  If  only  non-spread  ionograms  are  considered,  the  statistics 
of  all  months  improve  to  96%. 

A  clarification  about  our  definition  of  MUF(300Q)  is  in 
place  here.  In  the  BISA  program  the  F  region  h'(f)  curve  is 
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ERROR  DISTRIBUTION  OF  MUF (3000)  C (MANUAL  MUF  -  AUTO  MUF) /MANUAL  MUF3 

JANUARY  1980  GOOSE  BAY,  LABRADOR 


ERROR  DISTRIBUTION  OF  MUF(3000)  [(MANUAL  MUF  -  AUTO  MUF) /MANUAL  MUF] 

JULY  1980  GOOSE  BAY,  LABRADOR 


transformed  into  an  oblique  ionogram  by  multiplying  each  frequency 
with  the  transmission  factor  M(h') 

fob  ■  M<h')fverr 

The  transmission  function  M(h')  is  calculated  for  a  distance  of  3000 
km  by  fitting  a  polynomial  to  the  URSI  specified  data  set  [URSI  Hand¬ 
book  of  Ionogram  Interpretation  and  Deduction,  Second  Edition,  Nov. 
1972,  p.  21;  World  Data  Center  A  Report  UAG-23].  The  maximum  f  b  thus 
obtained  is  defined  as  MUF(3000)  in  the  present  BISA  program.  This 
may  be  either  the  MUF  FI  or  the  MUF  F2  whichever  is  larger.  For  the 
error  statistics  the  comparison  was  therefore  made  with  the  larger  of 
the  MUF  Fl/MUF  F2  values.  It  might  be  advisable  in  the  future  to 
autoscale  always  MUF  F2,  which  could  easily  be  do*e.  The 
significantly  higher  percentage  rates  for  the  non-spread  ior.ograms 
indicates  the  uncertainty  in  the  MUF  definition  for  spread  ionograms. 


5.1.3  M(3000) 

The  M(3000)  propagation  factor  is  derived  from  the 

MUF (3000) ,  as  defined  in  5.1.2,  by  dividing  it  by  foF2.  The  error 

distribution  functions  are  shown  in  Figures  12  to  15.  The  minimum 
accuracy  requirement  of  AM  *  ±0.2  is  fulfilled  for  82%  of  all 
ionograms.  In  January,  which  had  the  highest  magnetic  activity,  only 
78%  of  all  ionograms  pass  the  0.2  error  test.  For  the  other  three 
months  the  percentages  are  above  81%. 


5.1.4  Average  Diurnal  Variation  of  the  Scaling  Accuracy 

The  percentage  of  ionograms  that  were  scaled  with  an  foF2 
error  of  less  than  1  MHz  is  plotted  as  function  of  local  time  in 
Figures  16  to  19.  Not  surprisingly,  the  best  results  are  obtained  for 
the  daytime  ionograms  which  in  the  average  are  less  disturbed. 
January  shows  a  marked  difference  at  night  between  the  percentage 
curves  for  all  and  the  non-spread  ionograms.  The  same  behavior  is 
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displayed  by  the  MUF  error  curves  (Figures  20  to  23)  and  the  M-factor 
error  curves  (Figures  24  to  27). 


5.2  Other  Parameters  and  Statistical  Summary 

Table  1  summarizes  the  results  for  the  parameters  foFl , 
fminF,  fmin,  h'F  and  h'F2.  No  FI  layer  was  observed  in  January.  One 
can  see  that  only  69%  of  all  ionograms  are  scaled  within  an  h'F2  error 
of  10  km.  This  is  caused  by  the  fact  that  BISA  simply  lets  the  trace 
point  h ' ( f  m )  with  the  lowest  height  for  f  >  foFl  determine  h'F2.  The 
accuracy  of  this  parameter  can  be  improved  by  precisely  determining 
the  leading  edge  of  the  echo  trace  at  f  *  f  in  the  raw  ionogram. 
This  method  will  also  be  applied  to  the  h'F  height.  Table  2  lists  the 
results  for  the  E  region  parameters.  The  data  are  separated  into  day 
and  night  ionogram  to  better  see  the  effect  of  the  occurrence  of  night 
E  [King,  1962]  which  is  contained  in  the  foE  column.  Since  night  E 
occasionally  occurs  already  in  the  late  afternoon  the  day/night 
transition  was  made  when  the  predicted  foE  goes  below  2.5  MHz.  For 
more  than  92%  of  the  day  ionograms  the  foE  value  is  scaled  within  0.2 
MHz;  the  corresponding  value  is  71%  for  night  E.  For  an  error  limit 
of  0.5  MHz,  88%  of  the  night  ionograms  are  successfully  scaled.  The 
autoscaled  h'E  values  suffer  somewhat  because  the  E  trace  algorithm 
finds  foE  and  h'E  simultaneously,  optimizing  foE  at  the  account  of 
h'E.  This  parameter  will  be  improved  by  the  same  method  discussed 
above  for  h'F2.  The  critical  frequency  of  the  sporadic  E  layer  is 
autoscaled  within  0.5  MHz  of  the  manual  value  for  84%  of  the  day  and 
69%  of  the  night  ionogram.  The  auto  values  for  foEs  are  consistently 
lower  than  the  manual  readings.  The  reason  is  that  BISA  searches  for 
a  continuous  trace  of  vertical  0-echoes  and  terminates  the  trace  when 
there  are  no  0-echoes  within  300  kHz.  The  manual  scaler  often  extends 
the  trace  if  the  0-gaps  are  filled  with  oblique  and  X-echoes. 

Finally,  Table  3  is  a  summary  of  the  results  discussed 
earlier  for  foF2,  MUF(3000)  and  M(3000). 
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Four  Months  947  93.6  96.9  91.1  95.7  78.0  90.6 


6.0  HARDWARE  IMPLEMENTATION  -  ARTIST 


The  ARTIST  (Automatic  Real  Time  Monogram  Scaler  with  True 
height  calculation)  is  an  8086  microprocessor  based  data  cruncher 
which  is  contained  in  a  19"  rack  mountable  chassis.  The  required  high 
processing  speed  is  achieved  through  the  use  of  the  8087  arithmetic 
chip.  The  ARTIST  has  500K  of  random  access  memory  and  2  Mbyte  of 
floppy  disk  storage.  FORTRAN  77  is  used  as  source  language.  At  the 
present  time,  Digisonde  128PS  ionograms  are  read  from  magnetic  tape  (9 
tracks,  1600  bytes/inch)  and  processed  to  produce  the  traces, 
parameters  and  electron  density  profiles  as  discussed  in  the  previous 
sections.  Data  transfer  from  the  Digisonde  256  uses  DMA  techniques. 
To  interface  the  ARTIST  to  the  Digisonde  128PS  requires  some  hardware/ 
software  modifications.  Data  records  from  the  DGS  128PS  will  be 
accepted  by  the  ARTIST  via  an  1/0  port. 
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AUTOMATIC  CALCULATION  OF  ELECTRON  DENSITY  PROFILES 
FROM  DIGITAL  IONOGRAMS 
3.  PROCESSING  OF  BOTTOMSIDE  IONOGRAMS 
Bodo  W.  Reinisch  and  Huang  Xueqin* 
University  of  Lowell  Center  for  Atmospheric  Research 
Lowell,  Massachusetts  01854 

ABSTRACT 


Automatic  scaling  of  bottomside  Digisonde  ionograms 
gives  the  E  and  F  region  echo  traces  with  high  accuracy. 
Polarization  and  incidence  angle  information  in  the  ionogram 
enable  the  extraction  of  the  vertical,  ordinary  polarization 
echo  trace  from  quiet  as  well  as  disturbed  ionograms.  The 
scaling  algorithm  is  tested  with  Digisonde  ionograms  from 
Goose  Bay,  Labrador,  which  show  spread  F  about  50%  of  the 
time.  Despite  these  disturbed  conditions  foF2  is  determined 
within  one-half  megahertz  for  475  ionograms  out  of  577  during 
January  1980.  A  profile  inversion  algorithm  calculates  the 
electron  density  profile  from  the  autoscaled  h'(f)  points. 
Parabolic  profile  shapes  are  assumed  for  the  E-region  and 
the  valley  between  E  and  F-layer.  The  F-layer  is  approximated 
by  a  single  sum  of  Chebyshev  polynomials,  and  the  entire  pro¬ 
file  is  described  by  a  set  of  13  numerical  values.  The  CPU 
time  required  for  ionogram  scaling  and  electron  density  pro¬ 
file  calculation  is  14  seconds  on  a  Cyber  71  computer.  In 
addition  to  the  profile  the  program  outputs  all  relevant 
ionogram  parameters:  critical  frequencies,  minimum  virtual 
heights,  frequency  and  range  spread,  average  echo  amplitudes 
and  Doppler  shifts. 

*0n  leave  from  Research  Institute  of  Radiowave  Propagation, 
Xinxiang,  Henan,  China. 
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1.  INTRODUCTION 


Automatic  scaling  of  vertical  ionograms  from  the 
many  ground-based  ionosondes  has  long  been  a  desirable  goal. 
Only  nowadays  when  advanced  digital  sounders  like  the  Digi- 
sonde  [Bibl  and  Reinisch,  1978]  are  available  can  this  goal 
be  reached.  Development  of  adequate  sounding  techniques 
was  the  unconditional  prerequisite  for  the  formation  of  suc¬ 
cessful  algorithms  to  automatically  scale  the  E  and  F  region 
echo  traces  in  a  vertical  ionogram.  In  the  past,  some  very 
impressive  results  were  occasionally  shown  where  autoscaling 
methods  were  applied  to  a  few  selected  ionograms.  What  we 
want  to  show  in  this  paper  is  an  autoscaling  algorithm  that 
works  with  quiet  and  disturbed  ionograms,  and  that  vertical 
electron  density  profiles  can  routinely  be  determined  at 
each  modern  Digisonde  station.  Ionograms  from  the  Goose  Bay 
Ionospheric  Observatory  (64.6N  geomagnetic)  of  the  Air  Force 
Geophysics  Laboratory  form  the  data  base  of  the  study.  These 
ionograms  display  a  large  variety  of  features:  quiet  and 
disturbed  daytime  recordings,  heavy  spread  F  during  the  night, 
the  mid-latitude  trough  moving  over  the  station,  fast  varia¬ 
tion  of  the  ionospheric  parameters  and  frequent  absorption 
events . 

The  typical  features  of  a  quiet  daytime  ionogram 
are  shown  in  Figure  1:  a  weak  E-trace,  transition  to  the 
F-trace,  splitting  in  0  and  X  polarization  echoes  and  double 
echoes.  The  lower  part  of  the  figure  shows  the  amplitude 
ionogram,  the  upper  part  the  so-called  status  ionogram,  ex¬ 
plaining  the  echo  type  for  each  amplitude  pixel.  The  iono¬ 
gram,  distinguishes  0  and  X  echoes,  vertical  and  off-vertical 
echoes  as  well  as  positive  and  negative  Doppler  shifts.  The 
scaling  algorithm  extracts  the  vertical  ordinary  echo  trace 
from  which  the  electron  density  profile  is  calculated. 


Section  2  describes  the  algorithm  for  the  autoscal¬ 
ing  of  the  bottomside  Digisonde  ionograms,  and  Section  3 
explains  the  profile  inversion  technique  for  autoscaled  data. 
Reference  is  frequently  made  to  the  two  earlier  papers  in 
this  sequence  of  three  [Reinisch  and  Huang,  1982  and  Huang 
and  Reinisch,  1982],  which  are  cited  as  Part  1  and  Part  2. 

2.  AUTOSCALING  OF  DIGITAL  BOTTOMSIDE  IONOGRAMS 


It  has  long  been  recognized  that  a  network  of  bot¬ 
tomside  ionosondes  provides  a  wealth  of  ionospheric  informa¬ 
tion  valuable  for  research  and  communication.  The  tedious 
effort  involved  in  scaling  and  evaluating  the  conventional 
analog  ionograms  has  been  a  consistent  deterrent  to  wide¬ 
spread  use  of  ionogram  data.  In  recent  years  digital  iono¬ 
sondes  became  available  and  a  serious  effort  began,  to  find 
ways  to  automatically  scale  the  ionograms.  At  that  point 
it  seemed  tempting  to  no  longer  generate  complete  ionograms 
but  to  reduce  the  soundings  to  a  few  frequencies  centered 
around  the  predicted  critical  frequencies  of  the  E,  FI  and 
F2  region  and  to  merely  retrieve  a  few  critical  parameters. 
The  complexity  of  disturbed  ionograms  especially  in  the  auro¬ 
ral  and  equatorial  regions  makes  this  approach  all  but  impos¬ 
sible.  Even  a  complete  digital  ionogram  which  merely  records 
amplitudes  as  function  of  range  and  frequency  (lower  part 
of  Figure  1)  is  inadequate  for  autoscaling.  The  amplitude 
information  must  be  complemented  by  polarization  flagging, 
and  the  distinction  between  vertical  and  off-vertical  echoes. 
The  Digisonde  at  Goose  Bay  routinely  records  these  multi¬ 
parameter  ionograms  on  magnetic  tape.  This  made  it  possible 
to  develop  ionogram  processing  software  at  a  main  frame  com¬ 
puter  (Cyber  71  or  CDC  6600)  and  test  its  performance  by 
applying  the  software  to  a  large  number  of  ionograms. 
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Signal-to-Noise  Enhancement 

The  location  of  a  station  determines  the  noise 
level  in  the  ionograms ,  since  it  is  generally  man-made  inter¬ 
ference  rather  than  white  thermal  or  atmospheric  noise  which 
obscures  the  echo  traces.  Although  the  phase  coherent  spec¬ 
tral  integration  in  the  Digisonde,  which  transmits  high-fre¬ 
quency-phase  coded  pulse  signals,  significantly  reduces  un¬ 
wanted  signals,  strong  interferers  still  contaminate  the 
ionograms  in  form  of  discrete  noise  lines.  The  typical 
signature  of  interference  is  a  fairly  constant  distribution 
of  amplitudes  over  the  range  bins,  while  the  echo  signals 
occupy  only  two  to  three  range  bins.  Because  of  the  dis¬ 
crete  nature  of  the  man-made  interference  the  data  for  each 
ionogram  frequency  is  processed  individually.  The  most 
probable  amplitude  found  over  all  range  bins  for  a  given 
frequency  is  used  to  set  the  noise  threshold.  During  severe 
spread  F  conditions,  the  most  probable  amplitude  could  be 
the  signal  level,  setting  too  high  a  threshold.  To  avoid 
this  problem,  two  probability  curves  are  formed,  one  for 
the  lower  64  and  one  for  the  higher  64  range  bins.  (The 
Digisonde  128PS  uses  128  range  bins.)  The  distribution  that 
peaks  at  a  lower  amplitude  sets  the  threshold  [Reinisch  and 
Smith,  1976]. 


E-Pegion  Trace 

The  task  of  determining  the  h’(f)  function  for 
the  vertical  echoes  with  ordinary  polarization  is  divided 
into  the  two  parts  of  finding  the  E-trace  and  the  F-trace. 
Since  one  can  devise  a  generic  function  to  approximate  the 
normal  E-trace,  which  is  impossible  for  the  F-trace,  it  is 
a  relatively  easy  task  to  find  the  E-trace.  Assuming  a  para¬ 
bolic  electron  density  distribution  in  the  E-region  [Bradley 
and  Dudeney,  1973] 
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(1) 


N(z)  =  Nr  {1  -  C-rr-=02} 

L  yE 

one  can  calculate  the  virtual  height  function  h*(f): 
zCf ) 

h»(f)  =  J  y*  dz;  (2) 

0 

N(z)  =  electron  density  at  height  z 

Ne  =  electron  density  at  height  zE,  the  peak  of  the  E-layer 

yE  =  half  thickness  of  the  parabolic  profile 

f  =  transmitted  frequency 
y 1  =  group  index  of  refraction 

fE  =  foE  =  9.0  MKS  units,  critical  frequency  of  E-layer. 

Neglecting  the  geomagnetic  field  it  is  easy  to  analytically 
calculate  the  virtual  height  as  [Davies,  1965,  p.  135]: 

h'(f)  =  zE  -  yE  +  7  yE  .  (3) 

A  simple  search  algorithm  determines  the  three  unknown  param¬ 

eters  fE,  zE  and  yE>  Each  h*(f)  function  defined  by  equation 
(3)  traces  out  a  sequence  of  amplitude  pixels  in  the  iono- 
gram.  The  sum  along  the  trace  of  all  amplitudes  with  ordi¬ 
nary  and  vertical  identification  is  the  measure  for  the  qual¬ 
ity  of  the  fit.  The  parameters  fE,  zE  and  yE  are  varied 
until  the  maximum  sum  is  found.  Actually,  the  amplitude 
sum  is  always  normalized  by  dividing  it  by  the  number  of 
frequency  increments,  and  this  average  amplitude  is  optimized 
The  optimum  trace  follows  the  peaks  the  echo  pulses.  The 
entire  trace  is,  therefore,  lowered  until  the  average  ampli¬ 
tude  per  frequency  is  reduced  by  6  dB.  This  lowered  trace 
follows  the  leading  edge  of  the  echo  pulses,  and  can  directly 
be  compared  with  the  manually  scaled  h’Cf)  trace. 
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Fortunately,  the  range  of  the  parameter  values  is 
rather  small.  The  critical  frequency  f£  can  be  predicted 
within  ±0.3  MHz,  i.e.  only  seven  f^.  values  have  to  be  tested 
for  ionograms  with  100  kHz  frequency  increments.  The  half 
thickness  y^  is  varied  from  10  to  40  km  in  5  km  steps,  re¬ 
sulting  in  seven  values.  The  height  of  the  layer  peak,  z£, 
is  varied  from  95  to  150  km,  meaning  12  values  in  case  of 
5  km  range  increments.  The  total  number  of  h*(f)  functions 
to  be  tested  is7*7xl2=588.  The  status  identification 
of  the  echo  amplitude  pixels  allows  this  simple  fitting  pro¬ 
cedure,  otherwise  X-echoes  would  cause  severe  errors.  The 
method  works  very  reliably  even  under  disturbed  conditions. 

In  Figures  2  and  3  the  letters  0,  X  and  B  indicate  vertical 
ordinary  and  extraordinary  echoes  and  oblique  echoes,  re¬ 
spectively.  The  autoscaling  (dark  circles)  properly  follows 
the  0-echoes. 

Special  provisions  for  high  latitude  stations  allow 
the  detection  of  particle  ionization  in  the  E-region,  often 
called  night  E  [King,  1962].  As  soon  as  the  predicted  foE 
value  goes  below  2.5  MHz,  the  program  tests  for  critical 
frequencies  of  up  to  6  MHz  or  the  center  (see  below)  of  the 
F-trace,  whatever  is  lower  in  frequency.  Figure  3  is  an 
example  of  a  night  E  condition  in  Goose ‘ Bay  with  foE  =  2.6 
MHz  at  20:19  local  time.  Ordinary  echo  traces  beyond  foE 
are  identified  as  Es  (foEs  =  4.0  MHz  in  Figure  2  and  3.1 
MHz  in  Figure  3). 


F-Region  Trace 

Part  1  discussed  the  methods  of  autoscaling  the 
F-region  echoes  for  topside  ionograms,  and  very  similar  meth¬ 
ods  apply  to  the  bottomside  soundings.  The  large  variations 
in  height  and  shape  of  the  F-region  traces,  together  with 
the  occurrence  of  multiple  echoes  and  spread  F  make  it  impos¬ 
sible  to  use  a  simple  curve  fitting  procedure  for  the  entire 
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F- trace.  During  disturbed  ionospheric  conditions,  it  is 
often  difficult  to  identify  the  main  overhead  echo  trace. 

In  Goose  Bay  this  is  especially  true  when  the  mid-latitude 
F-region  trough  [Sharp,  1966;  Bowman,  1969]  moves  over  the 
station.  The  direction  finding  capability  of  the  Digisonde 
identifies  the  bulk  of  the  oblique  echoes,  but  angular  res¬ 
olution  and  sidelobes  of  the  antenna  array  pattern  can  cause 
erroneous  identifications.  In  analogy  to  visual  inspection, 
the  automatic  trace  identification  starts  by  finding  the 
center  of  the  F-trace,  i.e.  the  region  where  the  change  of 
h'  with  frequency  is  small  and  the  echo  amplitude  is  strong. 

A  first  rough  trace,  called  the  baseline,  is  then  constructed 
by  sliding  a  searching  window  (5  frequencies  x  35  range  bins) 
frequency-by-frequency  from  the  trace  center  toward  higher 
and  lower  frequencies.  In  an  effort  to  minimize  the  required 
CPU  time,  the  program  uses  selected  data  points,  called  first 
local  maxima  (FLM).  For  each  frequency  the  position  of  the 
first  (in  the  direction  of  increasing  range)  amplitude  maxi¬ 
mum  of  an  O-echo  between  two  zero  amplitudes  is  defined  as 
an  FLM.  The  number  of  FLM's  per  frequency  is  not  limited. 
Each  FLM  is  checked  whether  it  is  a  double  or  triple  echo 
by  comparing  its  amplitude  with  the  pixel  at  half  or  two- 
thirds  the  range.  If  the  FLM  is  smaller  in  amplitude,  it 
is  deleted. 

The  frequency-by-frequency  searching  process  de¬ 
fines  a  rough  trace  close  to  the  leading  edge  of  the  pulse 
and  is  able  to  successfully  trace  the  F1/F2  transitions, 
but  it  also  produces  some  undesired  h*  discontinuities  when 
man-made  interferers  distort  the  echo  pulses  at  certain  fre¬ 
quencies.  Also,  in  the  F2  cusp  region,  the  baseline  is  gen¬ 
erally  a  poor  representation  (Figure  3)  of  the  ordinary  F- 
trace  when  spread  F  and  external  interference  are  present. 

Since  the  ordinary  vertical  echoes  are  identified 
in  the  Digisonde  ionograms,  piecewise  smoothing  of  the  data 
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is  possible,  in  contrast  to  the  situation  described  in 
Part  1  for  the  topside  ionograms.  A  six-point  smoothing 
method  based  on  the  signal  amplitudes  is  described  in  detail 
in  an  unpublished  report  [Reinisch  et  al,  1982],  For  fre¬ 
quencies  below  the  cusp  region  the  smoothed  baseline  runs 
just  below  the  peaks  of  the  echo  pulses.  The  fast  rise 
in  h1  at  the  F2  cusp  and  the  frequent  occurrence  of  spread 
F  require  a  special  fitting  procedure  for  the  trace  near 
foF2 . 


Hyperbolic  Trace  Fitting  at  the  F2  Cusp 

A  polynomial  or  a  suitable  rational  function 
can  describe  the  F2  cusp  with  reasonable  accuracy.  The 
simple  hyperbolic  function 

hi  =  r  *  alET  a  >  0,  b  <  0  (4) 

has  a  pole  at  f  =  -  a/b  and  is,  therefore,  well  suited 
to  fit  the  trace  near  the  critical  frequency.  The  coeffi¬ 
cients  r,  a  and  b  are  determined  in  such  a  way  that  the 
h' -function  follows  the  pixels  with  the  maximum  amplitudes, 
assuring  that  a  reasonable  cusp  will  be  defined  even  in 
the  presence  of  spread  F.  Frequently  the  0-trace  near 
foF2  is  not  visible  on  the  ionograms  because  of  ionospheric 
conditions  or  high  interference,  while  the  X-trace  is 
clearly  developed.  The  reliability  of  the  autoscaling 
is  therefore  largely  improved  if  an  0  and  an  X  hyperbola 
are  simultaneously  fitted  to  the  polarization  tagged  ampli¬ 
tude  data.  If  equation  (4)  is  the  X-hyperbola,  the  0- 
hyperbola  is  given  by 

K  -  r  *  a+EF*"  <5) 

with  f*  =  f  +  ^  fH»  where  the  gyrofrequency  fH  is  evaluated 
for  300  km  altitude.  The  twc  hyperbolas  (4)  and  (5)  are 
simply  offset  by  half  the  gyrofrequency,  and  they  are  posi- 


tioned  on  the  ionograms  such  as  to  maximize  the  sum  of  the 
0  and  X  amplitudes  in  narrow  stripes  around  the  0  and  X- 
hyperbolas,  as  shown  in  Figure  4.  For  each  hyperbola  point 
the  amplitudes  are  summed  downward  to  the  height  of  the 
preceding  frequency.  This  method  forces  the  hyperbola 
to  the  inner  edge  of  the  trace  cup.  The  computer  printout 
of  the  ionogram  in  Figure  4  does  not  show  the  signal  ampli¬ 
tudes,  but  instead  the  polarization  tags  0  and  X  for  the 
vertical  echoes,  and  the  letter  B  for  the  oblique  echoes. 
Only  O-signals  contribute  to  the  0-hyperbola’ s  amplitude 
summation,  and  X-signals  to  the  X-hyperbola's.  The  hyper¬ 
bolic  traces  for  the  ionogram  in  Figure  4  are  marked  by 
heavy  dots  (O-trace)  and  X’s  (X-trace).  The  hyperbolic 
section  of  the  trace  ends  where  the  O-hyperbcla  intersects 
the  baseline  (3.3  MHz),  and  the  baseline  then  defines  the 
trace  toward  lower  frequencies. 

An  adequate  set  of  hyperbolas  in  the  neighborhood 
of  the  cusp  region  must  be  tried  to  find  one  that  simulates 
the  F  cusp.  Using  the  baseline  points  as  a  guide  the  method 
of  min-max  rational  functions  CScheid,  1968,  p.  289]  fits 
hyperbolas  to  subsets  of  three  points.  Starting  at  the 
high  frequency  end  of  the  baseline,  a  hyperbola  is  fitted 
to  the  last  three  baseline  points  by  minimizing  the  distance 
x  from  the  hyperbola  (see  Figure  5).  The  range  deviations 


of  the  three  points  (f..,  ht)  from  the  curve  are 

h!  -  {r  +  . }  =  x  i  =  1,  2,  3.  (6) 

These  equations  can  be  written  in  the  form  (deleting  the 
primes  on  the  heights): 

(hj^  -  r  -  x)a  +  (h^  -  r  -  x)f^  b  -  1  =  0  (7) 

( h 2  -  r  +  x)a  +  (h2  -  r  +  x)f2  b  -  1  =  0  (8) 

(h3  -  r  -  x)a  +  (hj  -  r  -  x)fg  b  -  1  =  0  (9) 
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Eliminating  a  and  b  leads  to  a  quadratic  equation  for  the 
determination  of  x 

-  r  -  x  (h^ 

h2  -  r  +  x  (h2 

hg  -  r  -  x  (hg 

or, 

X  =  -  Ijr  (1  -  /I  -  —  (11) 

ZA 


r  -  x)  f1  1 

r  +  x)  f2  1 


r  -  x)  fg  1 


(10) 


where 

A  =  2  <fx  -  f3) 

B  =  Bq  +  Ar 

C  =  C0  "  C1  r 

Bq  =  -  f1(h3+2h1-h2)  -  f2^h3-h1)  -  f3(h2-2h3-h1) 
Cq  =  f1h1(h3-h2)  +  f2h2(h1-h3)  •>  f3h3(h2-h1) 

0.^=  f1(hg-h2)  +  f2(h^-h^J  +  f3(h2-h1). 


(12) 

(13) 

(14) 

(15) 

(16) 
(17) 


The  coefficients  B  and  C  are  functions  of  the  base  height  r, 
which  is  unknown.  Inspection  of  equation  (11)  shows  that  |x| 
will  be  a  minimum  if  4AC/B2  =  0.  This  may  not  be  possible 
to  achieve,  but  one  can  find  the  minimum  value  of  | AC/ B 2 | , 
by  setting  the  derivative  with  respect  to  r  equal  to  zero: 

(~)  =  0  (18) 
B 


resulting  in 


r 


B0  + 


2  AC, 


(19) 
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Once  r  and  x  are  calculated  from  (19)  and  (11)  the  coeffi¬ 
cients  a  and  b,  specifying  the  X-hyperbola,  can  be  found. 

The  corresponding  O-hyperbola  is  determined  by  simply  replac¬ 
ing  f  by  f*.  Finally  the  0  and  X  sums  are  determined. 

It  is  of  course  unlikely  that  the  hyperbolas  fitted 
to  the  last  three  points  (highest  frequency)  of  the  baseline 
will  be  a  good  approximation  for  the  traces  in  the  cusp 
region.  The  process  is  repeated  for  the  next  three  points 
and  so  on,  down  to  frequencies  at  the  center  window.  The 
same  process  is  then  applied  to  triple  points  spaced  by 
two  frequencies  and  by  three  frequencies.  By  numbering 
the  baseline  data  points  from  1  to  N,  one  can  write  the 
triple  point  sets  to  which  hyperbolas  are  fitted  in  the 


following  way: 

GROUP  1 

P1  P2 

P3 

GROUP  2 

P1  P3 

PS 

GROUP 

P1  P4 

3_ 

P7 

P2  P3 

P4 

P2  P4 

P6 

P2 

P5 

P8 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

PN-2  PN-1 

PN 

PN-4  PN-2 

PN 

PN-6 

PN- 

3  PN 

The  total  number 

of  trial  hyperbolas 

is  3N- 

12. 

From 

these 

hyperbolas  the  one  with 

the  maximum 

0  and  X 

amplitude 

sum 

is  selected  for  the  trace  representation.  Combining  the 
O-hyperbola  with  the  smoothed  baseline  results  in  a  prelim¬ 
inary  trace  for  the  O-polarization  echoes.  Lowering  the 
entire  curve  to  the  leading  edge  of  the  echo  pulses  gives 
the  final  F-trace.  Figures  6  and  7  illustrate  the  perform¬ 
ance  of  the  scaling  algorithm  for  quiet  and  disturbed  condi¬ 
tions.  Figure  7  shows  the  F-trace  before  it  had  been 
lowered. 
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Parameter  Extraction 


At  the  present  time,  the  program  determines  the 
following  ionospheric  parameters:  foF2,  foFl,  fminF,  MUF 
(3000),  M(3000),  range  spread  F,  frequency  spread  F,  h*F, 
foE,  foEs,  range  spread  E,  frequency  spread  E,  h*E,  h'Es 
(parameters  are  defined  in  the  U.R.S.I.  Handbook  of  Ionogram 
Interpretation  and  Reduction,  Second  Edition,  November 
1972,  by  W.  R.  Piggott  and  K.  Rawer,  World  Data  Center 
A,  Report  UAG-23).  In  addition,  averaged  echo  amplitude 
values  are  given  for  each  megahertz,  normalized  to  a  reflec¬ 
tion  height  of  100  km.  The  Doppler  signatures  of  the  echoes 
are  also  saved  for  further  analysis.  Statistical  evaluation 
of  these  autoscaled  parameters  will  be  discussed  in  another 
paper,  but  to  illustrate  the  performance  level  of  the  scaling 
algorithm  the  comparison  between  manual  and  autoscaled 
foF2  values  is  given  in  Figure  8.  The  dashed  curves  are 
the  error  distribution  for  577  hourly  ionograms  from  Goose 
Bay,  January  1980;  82%  of  these  ionograms  have  an  error 
of  less  than  0.5  MHz.  The  solid  curves  include  only  the 
256  ionograms  without  spread  F.  The  difference  between 
the  solid  and  the  dashed  curve  is  acceptably  small,  showing 
that  the  algorithm  performs  reasonably  well  also  in  the 
presence  of  spread  F.  as  demonstrated  in  Figure  7.  The 
autoscaling  of  the  F-trace  during  this  trough  affected 
ionograms  is  remarkably  accurate.  The  ionograms  show  several 
cusps  and  the  one  selected  by  the  program  is  the  vertical 
echo  trace  with  the  strongest  amplitudes.  Differences 
in  interpretation  in  this  type  ionogram  are  the  main  reason 
for  foF2  errors  of  1  to  2  MHz. 


3.  CALCULATING  THE  VERTICAL  ELECTRON  DENSITY  PROFILE 


Part  2  discussed  an  efficient  method  of  inverting 
the  hf(f)  trace  into  the  vertical  electron  density  profile. 


using  the  autoscaled  h*(f)  data  points  as  input.  The  profile 
is  described  by  a  sum  of  Chebyshev  polynomials  the  coeffi¬ 
cients  of  which  are  optimized  in  a  least-squares  sense  to 
best  reproduce  the  measured  h’(f)  trace.  This  profile-fitting 
method  is  well  suited  for  autoscaled  data  and  describes  a 
monotonic  profile  with  a  single  polynomial.  The  bottomside 
profile  may  have  an  ionization  valley  between  E  and  F  region 
and  the  polynomial  ansatz  must  be  expanded.  Figure  9  shows 
a  reasonable  model  of  a  vertical  electron  density  distribution, 
which  is  convenient  for  calculational  purposes:  up  to  the 
height  z the  profile  is  described  by  the  E-region  parabola 
defined  by  f^.,  y^.  and  z^,,  the  F-region  above  z2  is  represented 
by  a  single  polynomial  (as  in  the  case  of  the  topside  pro¬ 
file),  and  the  valley  around  zv  is  represented  by  another 
parabola  which  smoothly  connects  to  the  E-parabola  at  z^  and 
the  F-profile  at  z2.  While  this  valley  may  not  precisely 
describe  the  actual  profile  in  this  height  region  it  allows 
for  a  better  profile  presentation  than  the  assumption  of  a 
monotonic  distribution.  A  reasonable  D  region  profile  could 
be  added,  modifying  the  E-parabola  for  f  0. 


Model  Parameters 

The  E-region  parabola  was  already  determined  in  the 
scaling  algorithm,  so  f^,,  z^  and  y £  are  known.  If  the  pro¬ 
file  between  z,  and  z0  were  known  one  could  calculate  the 

X  L 

F-profile  with  the  profile-fitting  technique,  since  z2  =  zE  +  w 
could  be  derived.  Using  the  expressions  derived  in  Part  2 
one  can  write  for  the  contribution  of  the  F-region  to  the 
vertical  height  at  frequency 


>  f 


E* 


AhF(fk)  =  JQ  Ai  Sik‘ 


I  <  7, 


(20) 
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(21) 


f2  -  fN2 

t2  =  — -  or  fN2  =  f2  -  t2  (f2  -  f2). 

f  2  x;  2  K  K  L 

rk  “  rE 

At  the  F-layer  peak  g  =  0  and 
ZF  =  AI+1; 

ft 

at  height  z2  the  g  value  is  1,  and  T^(l)  =  1,  so 

I 

z2  =  AI+1  +  i|Q  Ai» 
or 


(24) 


(25) 


(26) 


I 

ZF  -z9=~  I  A  •  •  (27) 

t  1  i=0  1 

All  these  expressions  were  derived  in  Part  2  for  the  topside 
profile;  the  plasma  frequency  at  the  satellite,  fNg,  is  re¬ 
placed  for  the  bottomside  by  f^,.  Only  the  ordinary  echo 
trace  is  considered,  as  is  usually  done  for  the  bottomside 
ionosphere.  Dependence  on  the  height  variation  of  the 
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geomagnetic  field  is  less  critical  and  it  suffices  to  calcu¬ 
late  the  S^k  for  constant  gyrofrequency  fH  and  dip  angle  0 
at  a  height  of  300  km.  For  given  fE  and  fp  the  S^k  can  be 
calculated  for  all  scaled  frequencies  fk»  The  group  refrac¬ 
tive  index  y*  for  the  ordinary  polarization  was  given  in 
the  Appendix  of  Part  2. 

The  virtual  height  of  an  F-echo  at  frequency  fk 
is  composed  of  E,  valley  and  F-region  contributions  (see 
Appendix) : 

h*  =  f  t  =  Ah’E  +  Ah' y  +  Ah'p  (28) 

(c  =  speed  of  light  in  free  space,  t  =  observed  group  travel 
time). 

The  F  profile  is  determined  by  minimizing  the 
square  error  of  the  differences  between  the  observed  virtual 
F  heights  h'k  and  the  values  obtained  from  (28): 

e  =  ^  Ch'k  -  CAh’Ek  Sik)32.  (29) 

h’k  is  measured  and  Ah’Ek  can  be  calculated  (see  equation 
A2 ) .  The  square  sum  is  taken  only  over  the  data  points  of 
the  lower  F-region  (K'  is  the  frequency  index  that  is  asso¬ 
ciated  with  the  minimum  virtual  height).  This  part  of 
the  F-trace  is  more  sensitive  to  the  valley  shape  than  the 
trace  at  higher  frequencies.  It  is  shown  in  the  Appendix 
that  the  valley  contribution  depends  on  the  valley  width  w 
and  the  slope  B  at  Z2.  In  an  iterative  procedure  w,  B,  and 
the  A^  coefficients  are  determined  such  as  to  minimize  the 
error  e  in  (29)  (see  the  Appendix  and  also  Part  2  for  the 
calculation  of  A^).  The  requirement  of  smooth  transitions 
between  the  three  profile  functions  gives  the  valley  param¬ 
eters  z^,  zv,  yv,  z 2  =  z£  +  w  and  fy.  The  F-profile  is  spec¬ 
ified  by  the  nine  (or  less)  coefficients  Aq,  A-^,  ...» 

Aj+1.  The  order  I  of  the  polynomial  expansion  (21)  is  auto- 


matically  reduced  from  its  maximum  value  of  7  if  the  profile 
shape  allows  it  (Part  2).  For  archiving  or  remote  data 
transfer  one  should  allocate  space  for  sixteen  numerical 
values  which  give  a  detailed  description  of  the  entire  elec¬ 
tron  density  profile:  f£,  z£,  yE,  fy,  zv,  yy,  fp,  AQ ,  A1# 

*  *  *»  AI+1. 


Discussion  of  Electron  Density  Profiles 

It  should  be  noted  that  the  calculations  for  the 
E  and  valley  regions  neglected  the  geomagnetic  field.  In 
principle,  this  would  not  be  necessary,  but  it  simplifies 
the  calculations  without  introducing  significant  errors, 
considering  the  limited  accuracy  of  the  autoscaling. 

The  automatic  vertical  electron  density  profiles 
for  the  ionograms  in  Figures  2,  3  and  6  are  shown  in  these 
figures.  As  a  measure  of  the  quality  of  the  automatic  pro¬ 
file  technique,  the  auto  profiles  are  compared  with  the  cor¬ 
responding  "manual"  profiles,  obtained  by  manually  scaling 
the  ionogram  and  applying  the  conventional  lamination  method 
for  the  true  height  inversion.  These  comparisons  are  shown 
in  Figures  10,  11,  and  12.  The  daytime  profiles  in  Figure 
10  show  variances  of  about  10  km,  but  in  general  the  agree¬ 
ment  between  auto  and  manual  profile  is  very  good.  The  third 
curve  is  calculated  by  applying  the  lamination  inversion 
method  to  the  autoscaled  data.  This  curve  is  very  close  to 
the  auto  profile  indicating  the  validity  of  the  profile-fit¬ 
ting  method.  Similar  good  agreement  is  shown  by  the  night 
profiles  in  Figure  11  for  the  disturbed  ionogram  in  Figure 
3.  The  day  profile  in  Figure  12  has  an  ionization  valley 
between  the  E  and  F  layer,  which  of  course  cannot  be  deter¬ 
mined  by  the  standard  lamination  technique.  The  valley 
minimum  is  located  at  f v  =  2.1  MHz  and  zv  =  139  km;  the 
width  of  the  valley  is  36  km. 
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For  comparison  with  the  International  Reference 
Ionosphere  IRI  [Rawer  et  al,  1978  and  1981]  it  might  be 
desirable  to  express  the  profile  as  a  sum  of  Epstein-step- 
f unctions  [Booker,  1977].  This,  however,  need  not  be  done 
at  the  ionosonde  site.  The.  automatic  profile  technique 
described  in  this  paper  is  tailored  for  real  time  application 
in  a  microprocessor  based  Digisonde  chassis  which  will  be 
field  tested  within  the  next  few  months. 


4.  CONCLUSION 


Ground  based  digital  ionograms  with  polarization 
and  incidence  angle  identifications  can  be  successfully  auto- 
scaled  for  quiet  and  for  disturbed  ionospheric  conditions, 
and  the  corresponding  vertical  electron  density  profiles  can 
be  calculated.  By  installing  automatic  real  time  ionogram 
scaler  with  all  modern  digital  ionosondes  could  provide 
an  always  up-to-date  global  ionospheric  model,  assuming 
a  sufficient  number  of  these  stations. 
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APPENDIX 


Parabolic  E  and  valley  regions: 


z  <  z_  -  y 


2*2 

fN2  =  f2  ci  -  c— ^>2: 

**  yr 


ZE  "  y  —  z  —  Z1 


z-z„ 

fj  a  ♦ 
v  yv 


ZX  —  z  —  z2* 


Virtual  height  contributions  for 

1  fk  fk+fE 

Ah'E  =  z£  -  y£  +  j  yE  y~  An  ZE  ~  yE  -  z  -  ZE 

£  )c  E 


Ah’ vi  =  fk  /I  In 


Y+4B  D 


—  +  /f 2 

2  n  k 


f  2  +  DY 
E  Y+4B2D 


(Al) 


(A2  ) 


$  zE  ^  z  ^  Zl* 


For  z^  <  z  <  z2: 


(A3) 


Ah* V2  =  2fk  B  /5  [sin-1 


/fk  "  fE  +  D 


+  sm 


'(y+4B2D)(f2  -  f2  +  D) 


(A4  ) 


w  =  z„  - 


0-1  -  dfNli  =  d  f2  r,  + 

B  ‘^,z=z2  ^fv[1+(yv)] 


fv 

2  T  (z-zv)’ 

yv 


(A5 ) 
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where  w  is  a  measure  of  the  valley  width  and  B  defines  the 
slope  at  z2  where  the  valley  parabola  joins  the  F-profile. 
This  means  that  B  must  also  equal  the  slope  of  the  F-profile 
at  z2: 


B  = 


4  f|  In  (f£/fF 


-  I  A.  Cl  +  2  (2-  T.)  ]. 

)  i=0  1  ^  1  g=l 


(A6  ) 


The  total  contribution  of  the  valley  region  z^  <  z  £  z2  is 

Ah'y  =  Ah’y-^  +  Ahfy2  *  Ah ' y  fj»»  y^>  w,  B)  (A7) 


where  w  and  B  are  the  only  two  unknowns.  Setting  w  =  0 
and  B  =  0  one  can  calculate  a  set  of  initial  coefficients 
A^  using  the  least-squares  error  method: 


E  =  Ji  th'k  "  Uh'Ek  +  Ah’vk  +  Jo  Ai  Sik>]2  <A8) 

where  X  is  the  number  of  F-trace  data  points.  Minimizing 
this  error  means  solving  the  set  of  1+1  equations  for  the 
A, : 


k|1  (h'k  “  Ah'Ek  "  Ah’vk}  Sjk  =  Ai  Qij »  ^  =  °»  ni»  I(A9) 
with 

K 

Qij  =  l  Sik  Sjk  i,  j  =  0,  1,  111,  I.  (A10) 

k-1 

At  frequency  f^,  g  =  1  and  z  =  z2  and  the  equation  for  z 


z 


AI+1  + 


reduces  to 


(All) 
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z2  =  AI+1  +  i|Q  Ai 


since  T^(l)  =  1. 


This  determines  the  last  coefficient: 


AI+1  =  ZE  *  w  -  J0  Ai 


(A12  ) 


The  slope  B  is  now  calculated  from  (A6)  and  a  new  set  of  co¬ 
efficients  A^  is  determined.  In  this  iteration  loop  for  B 
only  the  part  of  the  F-trace  between  f^  and  f^,  is  used, 
where  h'F  has  its  minimum  at  f^, .  The  process  is  then  re¬ 
peated  for  w  *  1  km,  2  km,  etc.  The  parameter  pair  w,  B  which 
has  the  smallest  least-square  error  in  equation  (A8)  is  used 
to  describe  the  valley  profile  and  to  calculate  the  final  set 
of  coefficients  A^  for  the  F-profile. 

The  valley  parameters  z^,  Zy,  fy,  yv  and  z2  are  ob¬ 
tained  from  the  smoothness  requirements  for  the  profile  func¬ 
tion  at  z1  and  z^.  At  z  =  z2*. 


Z,-Zv  2 

f|  s  fj  Cl  +  (‘T-Z>  3 


(A13  ) 


dCfN2)  z  =  z. 


2CZ2  -  Zy)  f2 


(A14  ) 


This  leads  to: 


yv  =  2  B  fy  /F 


(A15 ) 


z2  =  Zy  +  2BD. 

The  smoothness  condition  at  z 
equation 


(A16 ) 


=  z^  requires  that  the  quadratic 


i 


(A17  ) 


Z  —  Z-p  2  .  Z“ZV  2 

fl  C1  -  (^T>  3  *  'v  tl  *  <-yT>  3 


has  only  one  solution  for  z.  This  gives 


Z1  = 


yv  * 

ZV  Y  +  ZE  C4) 

Y  +  (JL) 

rv 


(A18  ) 


z v  =  zE  +  /D[Y  +  (yv/fv)2]. 


(A19) 


Substituting  (A15): 


Zy  Y  +  HZg  B2D 
Y  +  4  B2D 


CA20) 


zv  =  z£  +  /D[Y  +  4B2D] . 


The  parameter  D  is  a  function  of  w  and  B;  from  (A16); 


(A21) 


w  3  z0  -  zF  =  zv  -  zv  +  2  BD  =  t''D[Y+4  B2D]  +  2BD  (A22) 


D  =  Y  V  4WB*  (A23) 

The  definition  for  the  virtual  height  contributions  for  the 
valley  is  (neglecting  the  geomagnetic  field): 


Ah'  =  /  y '  dz  *  / 


AT 


Ah’  *  / 


(A24 ) 


Sp  J  *  Z— —  Z  2 

b  A-  CjA)  [i-c— ^ — )  ] 

•>  y  p 
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Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Digisonde  Ionogram 

Real  time  Digisonde  128PS  ionogram  with  optic¬ 
ally  weighted  font  displaying  signal  amplitudes 
(bottom)  in  4  dB  increments  and  status  (top). 
Sixteen  level  status  indicates  3  incidence 
angles,  0  and  X  polarization  and  4  Doppler 
frequencies. 

Autoscaling  for  Relatively  Quiet  Ionogram 

Normal  E,  sporadic  E,  FI  and  F2  traces  are  ac¬ 
curately  determined.  The  autoprofile  shows 
only  minute  inflection  at  foFl  =  5.7  MHz. 

Night  E  and  Spread  F 

E  and  Es  trace  as  well  as  the  maximum  F-trace 
are  correctly  identified.  The  baseline, 
labeled  S,  is  shown  for  reference.  The  solid 
line  is  the  auto-profile.  The  F  region  fre¬ 
quency  spread  is  fxl  -  fxF2. 

Amplitude  Windows  for  Hyperbolas 

The  0  (or  X)  amplitudes  below  each  hyperbola 
point  down  to  range  of  previous  frequency  are 
summed;  for  last  frequency  range  window  ex¬ 
tends  to  upmost  range  bin. 

Hyperbolic  Fitting  with  Min-Max  Method 

Equal  distance  x  of  three  points  from  hyperbola 
is  minimized. 

Autoscaling  for  Quiet  Conditions 

Small  scaling  errors  occur  at  8.8  MHz  (0  and  X 
trace  crossover)  and  at  5.0  MHz  (reduced  trans¬ 
mitter  power).  Weak  E-trace  is  properly  iden¬ 
tified.  Auto-profile  shows  valley  between  E 
and  F. 
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Figure  7 


Figure  8 


Figure  9 


Figure  10 


Figure  11 


Figure  12 


Autoscaling  During  Trough  Conditions 

F-trace  (before  lowering)  identifies  the  strong 
est  vertical  echo  sequence  while  mid-latitude 
trough  moves  over  Goose  Bay. 

Error  Distribution  of  foF2 

Relative  deviation  of  auto  foF2  from  manually 
scaled  values  centers  around  zero.  Inclusion 
of  ionograms  with  spread  F  widens  the  distri¬ 
bution  function  only  slightly. 

Electron  Density  Profile  Model 

A  parabolic  valley  smoothly  joins  the  E-region 
parabola  at  z^  and  the  F-region  profile  at  z2- 

Daytime  Profile  for  Quiet  Conditions 

Goose  Bay  21  April  1980,  10:00  AST.  Compari¬ 
son  of  manual  and  automatic  methods. 

Nighttime  Profile  During  Disturbed  Conditions 

Goose  Bay  5  September  1980,  20:19  AST. 

Daytime  Profile  with  Ionization  Valley 

Goose  Bay,  9  January  1980,  10:00  AST. 

Lamination  method  cannot  show  the  valley. 
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Figure  2 
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Figure  4 
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AUTOSCALING  DURING  TROUGH  CONDITIONS 

GOOSE  BAY  7  JAN  1980 


ULCAR 
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ERROR  DISTRIBUTION  OF  FOF2  (MANUAL  FOF2-  AUTO  FOF2) 
GOOSE  BAY  JANUARY  1980 
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ELECTRON  DENSITY  PROFILE  MODEL 


PLASMA  FREQUENCY  (MHZ) 

DAYTIME  PROFILE  for  QUIET  CONDITIONS 
GOOSE  BAY  21  APR.  1980,  10:00  AST. 
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PLASMA  FREQUENCY  (MHZ) 

NIGHTTIME  PROFILE  during  DISTURBED  CONDITIONS 
GOOSE  BAY  5  SEP.  1980,  20:19  AST. 
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PLASMA  FREQUENCY  (MHZ) 

DAYTIME  PROFILE  with  IONIZATION  VALLEY 
GOOSE  BAY  9  JAN  1980,  10:00  AST. 


